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Abstract 

Various 2-particle direct photon-hadron correlation yields in p + p and Au + Au 
collisions at ^Jsnn = 200 GeV are presented. The per-trigger yield of direct photon 
CN ■ hadron pairs from direct-photon-jet correlations is obtained by a statistical subtraction 

of the decay photon pairs from inclusive photon-hadron sample. The decay photon per- 
| trigger yields are estimated from the measured 7r°-hadron by means of a Monte Carlo 

based calculation which takes into account decay kinematics and detector response. 
t— I ■ Under the assumption that the suppression is nearly px independent using a specific 

averaging scheme, we find a ratio of Au + Au to p + p per-trigger photon yields Iaa 
^ ' averaged over all available pt bins consistent with the single particle suppression level 

^ • Raa, which can be interpreted as a qualitative confirmation of the basic geometrical 

picture of jet suppression at RHIC. The application of the event by event photon 
isolation cuts in p + p results our highest precision measurement yet, and allows for 
precision studies of the baseline fragmentation function D(z) and also a determination 
of the apparent intrinsic kx, or non-zero transverse momentum of the original collision 
partons. With a model dependent extraction method, the average \J < k T > at this -y/s 
in p + p is found to be approximately 2.5-3 GeV, consistent with analysis of di-hadron 
(di-jet) correlations [lj. Finally we present a unique direct measurement of prompt 
. photons from jet fragmentation. 

§s : 1 Introduction 

o ■ 

J> , At the Relativistic Heavy Ion Collider (RHIC), experimental results from the have 

established the formation of hot and dense matter of a fundamentally new nature in 
Au+Au collisions at y / J/v r /v=200 GeV [2]. One of the most important probes of this 
dense colored matter is energy loss by hard partons leading to a suppression of normal 
jet production of hard (E >~ 1 — 2 GeV) particles. As a complement to di-hadron 
correlations, which can directly access di-jet production and their structure [3], [I], 
direct photon-hadron correlations can be used to study photon-jet production in the 
medium without various biases and complications. This is due to the fact that since 
it lacks color charge, the photon escapes the dense QCD matter without interacting. 
For this reason, using high px prompt photons, which have been demonstrated to be 
unmodified when traversing the medium due to their lack of color charge [5] , as triggers 
in 7 - hadron (i.e. 7-jet) correlation studies has long been considered to be a "golden 
channel" to study energy loss. When the photon trigger is a product of the dominant 
Leading Order QCD Compton scattering process, the photon's 4-momentum should be 
the same as the opposite quark parton's. Therefore the photon in its primordial state 
can be used to directly measure the fragmentation functions of the opposite jet which 
itself is modified by final state effects in the case of Au+Au collisions. 
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The above picture is based on a Leading Order (LO) approximation and interpre- 
tation of events where there are always only two "clean" outgoing products, di-parton 
(g, q) or photon-parton for 7-jet, which are perfectly balanced in momentum in the 
direction transverse to the incoming colliding particles. In reality there are a num- 
ber of complications that could make this interpretation incorrect that need to be 
understood even in the baseline elementary p + p collisions. For example, Next to 
Leading Order (NLO) /fragmentation contributions of hard prompt photons in di-jet 
events obviously have different properties, in particular in Au + Au are expected to 
follow similar suppression effects as the di-jet events which makes making quantifying 
the baseline probability for their appearance in the overall 7-jet event samples crucial. 
Another example is the role of an apparent intrinsic net transverse momentum kx of 
the incoming colliding parton-parton system potentially due to initial state effects or 
non-perturbative or higher order gluon radiation, which at the value discovered in di- 
hadron correlations pQ can considerably change the actual momentum transfer (Q 2 ) 
being sampled in the kinematic ranges currently being studied. Thus, before studying 
7-jet production A+A, one needs to understand these intrinsic complications to the 
LO picture in p + p to gauge their effect on calculations and interpretations of energy 
loss observables. 

This analysis was performed using approximately 950 million Au+Au minimum 
bias events from the Run 4 data set and 471 million level- 1 photon-triggered events 
from the p+p Run 5 and Run 6 data sets. The Beam-Beam Counters (BBC) and 
Zero-Degree Calorimeters (ZDC) are used to trigger the minimum bias data. These 
detectors are also used to estimate the collision centrality. The p+p photon trigger 
requires that a module in the Electromagnetic Calorimeter (EMC) fire in coincidence 
with the BBC. 

The PHENIX central arms, each covering the 0.7 units of pseudo-rapidity around 
mid-rapidity and 90° in azimuth, contain charge particle tracking chambers and elec- 
tromagnetic calorimetry. The EMC consists of two types of detectors, six sectors of 
lead-scintillator sampling calorimeters and two of lead-glass Cherenkov calorimeters. 
The fine segmentation of the EMC (~ 0.01 x 0.01 for PbSc and ~ 0.008 x 0.008 for 
PbGl) allows for the reconstruction of the 27 ir° and r\ decays to pt > 20 GeV. 

Charged hadrons are detected using a drift chamber at a radial distance of 2.0 m 
and and a multi-wire proportional chambers (PCI) at a distance of 2.5 m. The momen- 
tum resolution was determined to be 0.7%® 1.0%p (GeV/c). Secondary tracks from 
decays and conversions are suppressed by matching tracks to hits in a second multi- 
wire chamber (PC3) and the EMC, both at distance of ~ 5.0 m. Track projections to 
the EMC plane are used to veto charged hadrons which shower in the EMC. 

Two-particle correlations are performed by measuring the azimuthal angle between 
photon triggers and charged hadrons. The correlation function 
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corrects for the limited acceptance of photon-hadrons pairs by dividing out the 
mixed events distribution. The correlation function is decomposed via the two-source 
model where the jet correlation is superimposed on an underlying event which is mod- 
ulated by elliptic flow. Hence, the jet function is expressed as J(A0) = C(Acp) — 
6o(l + 2(^2^2) cos ^ A0). The underlying event level, bo, is determined by the ZYAM 
procedure, described elsewhere [6]. 

We define a direct photon as any photon not from a decay process. It follows that 
the per-trigger yield (Y = l/N tr igger J{A<fi)) f° r direct photons may be obtained by a 
statistical subtraction of the decay per-trigger yield from the inclusive per-trigger yield 
according to: 
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where R-y= N inc i usive /N decay Ry is determined for Au+Au collisions in [5j and is 
derived from the 7r° [7J and direct photons spectra in p+p [8]. The direct photon 
yields from the statistical subtraction method do not, by definition, exclude photons 
from jet fragmentation or medium induced photon production. 

The decay photon per-trigger yields are determined from the parent (7r° or 77) 
per-trigger yields via a Monte Carlo procedure. A flat distribution of parent mesons 
are decayed into the PHENIX aperture uniformly in the z and phi directions. This 
determines the mapping of the parent to daughter pr, p(p^ ~^ Pt)i where p T is 
smeared by the detector resolution. In order to reproduce the correct prdependence 
of the decay photon distribution W is applied as a weight factor to the parent meson- 
hadron ^distribution on a pair-by-pair basis. The finite reconstruction efficiency of 
the parent mesons is corrected for using the published PHENIX px spectra [7J, [9]. 
The decay per-trigger yield from 7r°'s can then be expressed in terms of the parent 
per-trigger yield. 

2 Constraining Direct Photon-Jet Production 
in Elementary p + p Collisions 

2.1 Isolation Cut Analysis 

In addition to the above statistical subtraction method, standard photon isolation cuts 
were applied event by event in a new analysis in order to dramatically reduce the 
contamination of di-jet events with ir° decay or fragmentation photons in the 7-hadron 
event sample. To be considered isolated, the sum of p^from all tracks and and EMCal 
energies must be < 0.1^ in a cone around the photon of size AR = \J At/> 2 + Arj 2 = 
0.5. Statistical subtraction of the remaining contribution for isolated tt° production is 
achieved through a modified version of the statistical method above where isolated 7r° 
are used as input to the decay-photon calculation. 

As in the statistical method the analysis is performed as a function of Acf). We find 
good agreement with the statistical method results as shown in figure [TJ Since the sta- 
tistical subtraction method includes, in principle, the contribution from fragmentation 
photon triggers, the agreement between the two methods places a constraint on the 
magnitude of such a contribution. 



2.2 Per trigger yields and A wayside Jet Fragmentation 
Function Analysis 

Results from the isolation cut analysis are shown in figure [21 plotted in terms of the 
fragmentation variables zt = p9p soc /p7p and xe = zt cos A<f> (taking only the compo- 
nent of the associated hadron momentum that is in the same direction as the trigger 
photon). We measure in 6 photon trigger px bins covering the range 5-15 GeV/c. The 
inverse slope parameter of exponential fits to the xe distribution with 7r° and direct 
photon triggers are shown in figure [3l Comparisons with theoretical calculations are 
currently underway. 
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Figure 1: For two trigger photon bins as indicated on the figure, the awayside per-trigger 
yield vs. A<p for the statistical and isolation cut analysis. The two methods give consistent 
results with the isolation method having much improved precision due to having only a small 
decay photon background subtraction. 




Figure 2: Per trigger yields of direct photon-hadron pairs as a function of Xe/zt for various 
selections of trigger p?. The data have been scaled by factors of 10 for visibility. 
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Figure 3: The inverse slope parameter of exponential fits to the xe distributions for ir° direct 
photon triggers. 



2.3 Model Dependent Determination of kx 

Also in p], a method for extracting the apparent intrinsic kx itself was described and 
used to extract a value of kx= 2.68 ± 0.35 for 7r°-h correlations. Using isolated direct 
photon triggers, the method is simplified and can be also be used to extract a kx value 
in this channel. Please refer [TJ for details of the method. The simplification occurs 
due to the lack of an additional fragmentation variable on the nearside, i.e pxt = 
pxt = PT'y- The measured experimental variable p ou t whose distribution for various 
trigger px bins is shown in figured! is defined as px a sin A</>, or the associated hadron's 
transverse momentum component perpendicular to the trigger photon direction, and 
is proportional to the kx on an event by event basis. 
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Figure 4: distributions for various values of trigger px- The data have been scaled by 
factors of 10 for visibility. 

By finding the average p ou t for each trigger photon px bin as shown in [5] a) one 
can extract the quantity average kx/x where x is just the ratio of the true awayside 
jet momentum to that of the direct photon trigger. This can be extracted from other 
measurements as in [1] and eliminated, leaving a pure measurement of the kx- For now, 
we simply use the event generator PYTHIA 6.3 to generate direct photon processes 
with initial and final state radiation turned off, and phenomenological kx parameters 
as shown on [5] b) to extract x and make a comparison to PYTHIA, which shows that 
for PYTHIA- like distributions of x a kx parameter in the vicinity of 2.5-3 GeV gives 




Figure 5: a) y< p 2 out > vs trigger p T for direct photons b) k T /x vs trigger p T for direct 
photons along with calculations from PYTHIA 6.3. 



values similar to the data, well-consistent with the value of 2.68 found for di-hadron 
di-jet correlations. 

This implies that the same k? bias effects in the di-hadron correlations exist also in 
the direct photon-jet correlations. This makes the comparison with the perturbative 
QCD calculations, which at only NLO should not contain such kx bias modifications 
even more important, as well as understanding in energy loss models of Au + Au how 
such a large value of imbalance between the trigger photon and awayside jet might 
alter energy loss interpretations. 

2.4 Direct Measurement of Prompt Photons from Jet Frag- 
mentation 

PHENIX has made further strides towards understanding complications to the LO 
intepretations of direct photon-jet correlations, in studying the contribution of single 
prompt photons that occur in di-jet events, the so-called fragmentation photons pro- 
duced by NLO hard photon radiation or non-perturbative fragmentation processes. A 
direct measurement technique discussed in [TU] has been used to measure the angu- 
lar distribution of such photons with respect to trigger hadrons in events where high 
Pt photons are tagged to only be nearby such hard hadrons-an aniz-isolation selec- 
tion. Integrating this distribution one can find the fraction of such photons to total 
hadron-correlated photons from all sources including decay, shown in figure [6] which 
can possibly be related to the fraction of the total direct photon production rate due 
to these fragmentation photons. However due to the restricted kinematic region of the 
measurement, interpretations of this fraction and the angular distribution itself need 
input from sophisticated higher order pQCD calculations (which likely do not yet exist) 
and thus close attention from the theory community. Nonetheless, the measurement 
is exciting because it provides the first step towards making the measurement in the 
Au+Au environment where several interesting predictions of jet-medium enhancement 
of the rate of production of such photons exist. 

3 Constraining Energy Loss theory in Au+Au 
with 7- jet 

In previous conferences [lllfT2j PHENIX has presented Au+Au results of the statistical 
subtraction method for direct photon-hadron correlation yields. With the now statis- 
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Figure 6: Fraction of hadron-photon that contain a prompt (fragmentation) photon 
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Figure 7: Ratio Iaa (see text) for the lowest trigger direct photon p T bin. 



tically improved p + p results (for this section, using the pure statistical subtraction 
method in p+p for comparison, combining Run6 and Run6 statistics) and an expanded 
px range for the trigger photons in Au + Au, we can now divide the per-trigger yields 
for Au + Au and p + pin many pt bins making the ratio Iaa = ^^rect Indirect f° r an 
awayside integration range of 2ir/5 radians around ir. An example is shown for the 
direct photon (trigger) px bin 5-7 GeV/c in figure [7J It is apparent that there are very 
large uncertainties, but most points have positive yield in Au + Au and a value of Iaa 
between and 1 which indicates the basic expectation of suppression of the awayside 
jet. 

We find that different sources of uncertainties dominate in different ranges of trigger 
photon and associated hadron hadron px- As either px is increased the statistics 
in the raw correlation functions obviously decrease as the production probabilities 
follow steeply falling spectra. However the combinatoric background from pairs of hard 
particles and soft underlying or otherwise uncorrelated particles is reduced dramatically 
as the px of the associated hadron is increased, having the effect of reducing the 
systematic uncertainty from the underlying event subtraction. Further, as the px of 
the trigger photon is increased there is a larger fraction of direct photon-h pairs in the 
inclusive photon-h pair sample due to the increasing direct photon signal relative to 
the suppressed decay photon background, and thus the higher trigger photon px bins 
also have a reduction in the uncertainty from the subtraction in equation [TJ Because of 
these effects we find that the total uncertainty remains more constant with increasing 
px than the loss in statistical precision that the falling production probability would 
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Figure 8: a) Mean value Iaa for each trigger px bin b) over all px bins for 0-20%, 20-40%, 
and 40-92% centrality bins. 



normally manifest. For this reason, integrating over large px hadron or trigger px 
bins causes a loss of information since the steeply falling production probability causes 
uncertainties to be dominated by the source corresponding to the low end of the px 
bin. Therefore in order to effectively give the higher px bins more weight, we take the 
plain average of all Px,y>PT,h bins. We call this average the Mean Value Iaa akin to 
performing the functional mean value: 

MeanValue I AA = J I A A(Px SS ° C )dpT (2) 

Under the further assumption that Iaa remains constant with px this average 
indeed corresponds exactly to the functional mean value, and it is found that this 
assumption is satisfied to what is likely a sufficient degree (considering our large un- 
certainties which would dominate any such uncertainties in the assumption) both in 
measurements of di-hadron correlations [1] and [3] and in most theoretical predictions 
[131 [14] at sufficiently high px (associated px >— GeV/c). 

Results of the Mean Value Iaa are shown in figure The transition of different 
uncertainty sources from systematic to statistical is even more apparent in figure a) 
for central events, and indeed the Mean Value Iaa seems to favor a constant value 
vs. trigger direct photon px although the uncertainties are still much too large to rule 
out with any significance rather large possible trigger photon px dependence scenarios. 
Still, again taking the average of all Iaa points together, shown on figure [8] b), now 
also for more peripheral bins, we find that the Mean Value Iaa for the central events 
is significantly positive at a two-sigma level and consistent with the single particle 
suppression level Raa- 

This confirms the basic geometrical picture of energy loss at RHIC |15|.|16j because 
vertices in 7-h correlations, and the trigger direct photon themselves, come from the 
entire collision volume. Therefore the fraction of 7-/1 pairs that are observed without 
significant suppression, come from the same geometrical region that surviving single 
high px particles do. Dividing this quantity by the unsuppressed direct photon yield is 
analogous to dividing the surviving single particle production rate by the expectation 
from p + p multiplied over the entire production weighted geometric volume of the ini- 
tial distribution of all hard scattering production points, as is done in the construction 
of .R^This means that the away-side suppression in the 7 — jet channel should simply 
reflect the same geometry as the single particle picture and, if geometry plays a domi- 
nant role, give the same suppression level. Since this geometrical picture is believed to 
be a "surface bias" picture where the only di-jets initially produced near the surface 
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Figure 9: A0 distribution of direct photon-hadron correlations from Run7 data for 5-7 trigger 
direct photon bin. 



it could be a confirmation of the surface bias picture, although further comparison 
to the di-hadron Iaa comparing 7-hadron and di-hadron correlations together in the 
same energy loss framework, are necessary to make this statement, in order to rule out 
possible effects in the di-hadron Iaa that are believed not to be from geometry. 

Looking forward to the full release of our Run7 analysis, we also present a first look 
in figure [9] at the lowest trigger px bin using ~ 3.0 billion events, which entices us to 
look across several px bins to look for consistent behavior that may be consistent to 
jet-medium displaced peak as seen already in di-hadron correlations, although at this 
point uncertainties are too large to make a definitive statement with just a single pt 
bin combination. 



4 Conclusions 

Various 2-particle direct photon-hadron correlation yields in p+p and Au+Au collisions 
at ^/snn = 200 GeV have been presented. Under the assumption that the suppression 
is nearly px independent using a specific averaging scheme, we find a ratio of Au + Au 
to p + p per-trigger photon yields averaged over all available px bins, Iaa, consistent 
with the single particle suppression level Raa, which can be interpreted as a qualitative 
confirmation of the basic geometrical picture of jet suppression at RHIC. To the extent 
that the prevailing geometric picture is believed to be that of surface bias where only 
jets ejected near the surface dominately contribute to yields, [15], [16], it may be said 
to confirm this picture, however it should be noted that any geometric scenario would 
yield Iaa — Raa- Nonetheless future comparisons of Iaa from 7-jet to that of di-jets, 
may indeed yield some insight into the details of possible surface bias models. 

We have also presented the first event by event isolation cut 2-p correlation results 
at RHIC. The application of the event by event photon isolation cuts in p + p results 
our highest precision measurement yet, and allows for precision studies of the baseline 
fragmentation function D(z), and well as a variable p ou t which is proportional to the 
apparent intrinsic kx, or non-zero transverse momentum of the original collision par- 
tons. With a model dependent extraction method, the average \J < k T > at this y^s 
in p + p is found to be approximately 2.5-3 GeV, consistent with analysis of di-hadron 
(di-jet) correlations pQ. The possible implications of this and also the improved preci- 
sion in the isolated yields warrant detailed comparison with the baseline perturbative 
QCD (pQCD) calculations used in the various models of jet energy loss. Finally, we 
have presented a unique direct measurement of single prompt photons from jet frag- 



mentation, of both angular information and the fraction of these photons in the entire 
photon sample (including decay photons) in the vicinity of the jet cone and for specific 
Pt cuts. 
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Various 2-particle direct photon-hadron correlation yields in p + p and Au + Au 
ON I collisions at y / s/v r /v = 200 GeV are presented. The per-trigger yield of direct photon 

hadron pairs from direct-photon-jet correlations is obtained by a statistical subtraction 
^| , of the decay photon pairs from inclusive photon-hadron sample. The decay photon per- 

trigger yields are estimated from the measured 7r°-hadron by means of a Monte Carlo 
based calculation which takes into account decay kinematics and detector response. 
Under the assumption that the suppression is nearly pt independent using a specific 
averaging scheme, we find a ratio of Au + Au to p + p per-trigger photon yields Iaa 
averaged over all available px bins consistent with the single particle suppression level 
Raa, which can be interpreted as a qualitative confirmation of the basic geometrical 
D , picture of jet suppression at RHIC. The application of the event by event photon 

isolation cuts in p + p results our highest precision measurement yet, and allows for 
precision studies of the baseline fragmentation function D(z) and also a determination 
of the apparent intrinsic Jzt, or non-zero transverse momentum of the original collision 
partons. With a model dependent extraction method, the average \J < k^ > at this yfs 
in p + p is found to be approximately 2.5-3 GeV, consistent with analysis of di-hadron 
(di-jet) correlations [lj. Finally we present a unique direct measurement of prompt 

ON ■ photons from jet fragmentation. 
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1 Introduction 

On . 

, At the Relativistic Heavy Ion Collider (RHIC), experimental results from the have 

established the formation of hot and dense matter of a fundamentally new nature in 
Au+Au collisions at y / J/v r /v=200 GeV [2]. One of the most important probes of this 
dense colored matter is energy loss by hard partons leading to a suppression of normal 
jet production of hard (E >~ 1 — 2 GeV) particles. As a complement to di-hadron 
correlations, which can directly access di-jet production and their structure [3], [3J, 
direct photon-hadron correlations can be used to study photon-jet production in the 
medium without various biases and complications. This is due to the fact that since 
it lacks color charge, the photon escapes the dense QCD matter without interacting. 
For this reason, using high pr prompt photons, which have been demonstrated to be 
unmodified when traversing the medium due to their lack of color charge [5] , as triggers 
in 7 - hadron (i.e. 7-jet) correlation studies has long been considered to be a "golden 
channel" to study energy loss. When the photon trigger is a product of the dominant 
Leading Order QCD Compton scattering process, the photon's 4-momentum should be 
the same as the opposite quark parton's. Therefore the photon in its primordial state 
can be used to directly measure the fragmentation functions of the opposite jet which 
itself is modified by final state effects in the case of Au+Au collisions. 

The above picture is based on a Leading Order (LO) approximation and interpre- 
tation of events where there are always only two "clean" outgoing products, di-parton 



{g, q) or photon-parton for 7-jet, which are perfectly balanced in momentum in the 
direction transverse to the incoming colliding particles. In reality there are a num- 
ber of complications that could make this interpretation incorrect that need to be 
understood even in the baseline elementary p + p collisions. For example, Next to 
Leading Order (NLO) /fragmentation contributions of hard prompt photons in di-jet 
events obviously have different properties, in particular in Au + Au are expected to 
follow similar suppression effects as the di-jet events which makes making quantifying 
the baseline probability for their appearance in the overall 7-jet event samples crucial. 
Another example is the role of an apparent intrinsic net transverse momentum kx of 
the incoming colliding parton-parton system potentially due to initial state effects or 
non-perturbative or higher order gluon radiation, which at the value discovered in di- 
hadron correlations [1] can considerably change the actual momentum transfer (Q 2 ) 
being sampled in the kinematic ranges currently being studied. Thus, before studying 
7-jet production A+A, one needs to understand these intrinsic complications to the 
LO picture in p + p to gauge their effect on calculations and interpretations of energy 
loss observables. 

This analysis was performed using approximately 950 million Au+Au minimum 
bias events from the Run 4 data set and 471 million level- 1 photon-triggered events 
from the p+p Run 5 and Run 6 data sets. The Beam-Beam Counters (BBC) and 
Zero-Degree Calorimeters (ZDC) are used to trigger the minimum bias data. These 
detectors are also used to estimate the collision centrality. The p+p photon trigger 
requires that a module in the Electromagnetic Calorimeter (EMC) fire in coincidence 
with the BBC. 

The PHENIX central arms, each covering the 0.7 units of pseudo-rapidity around 
mid-rapidity and 90° in azimuth, contain charge particle tracking chambers and elec- 
tromagnetic calorimetry. The EMC consists of two types of detectors, six sectors of 
lead-scintillator sampling calorimeters and two of lead-glass Cherenkov calorimeters. 
The fine segmentation of the EMC (~ 0.01 x 0.01 for PbSc and ~ 0.008 x 0.008 for 
PbGl) allows for the reconstruction of the 27 7r° and rj decays to px > 20 GeV. 

Charged hadrons are detected using a drift chamber at a radial distance of 2.0 m 
and and a multi-wire proportional chambers (PCI) at a distance of 2.5 m. The momen- 
tum resolution was determined to be 0.7%® 1.0%p (GeV/c). Secondary tracks from 
decays and conversions are suppressed by matching tracks to hits in a second multi- 
wire chamber (PC3) and the EMC, both at distance of ~ 5.0 m. Track projections to 
the EMC plane are used to veto charged hadrons which shower in the EMC. 

Two-particle correlations are performed by measuring the azimuthal angle between 
photon triggers and charged hadrons. The correlation function 

C N mixed (A4>) 

corrects for the limited acceptance of photon-hadrons pairs by dividing out the 
mixed events distribution. The correlation function is decomposed via the two-source 
model where the jet correlation is superimposed on an underlying event which is mod- 
ulated by elliptic flow. Hence, the jet function is expressed as J(A0) = C(Acp) — 
6o(l + 2(^2 ^2 ) cos 2A0). The underlying event level, &o> is determined by the ZYAM 
procedure, described elsewhere [6]. 

We define a direct photon as any photon not from a decay process. It follows that 
the per-trigger yield (Y = l/N tr igger J(A0)) for direct photons may be obtained by a 
statistical subtraction of the decay per-trigger yield from the inclusive per-trigger yield 
according to: 

^ R-fYi nc l us i ve Ydecay /-. x 

" direct ^ i \ / 



where Ry= A T inclusive/^ T decay Rj is determined for Au+Au collisions in [5] and is 
derived from the 7r° [7J and direct photons spectra in p+p [8]. The direct photon 
yields from the statistical subtraction method do not, by definition, exclude photons 
from jet fragmentation or medium induced photon production. 

The decay photon per-trigger yields are determined from the parent (-7T or 77) 
per-trigger yields via a Monte Carlo procedure. A flat distribution of parent mesons 
are decayed into the PHENIX aperture uniformly in the z and phi directions. This 
determines the mapping of the parent to daughter px, p(Px ~^ Pt)i where p T is 
smeared by the detector resolution. In order to reproduce the correct prdependence 
of the decay photon distribution W is applied as a weight factor to the parent meson- 
hadron ^distribution on a pair-by-pair basis. The finite reconstruction efficiency of 
the parent mesons is corrected for using the published PHENIX px spectra [7], [9]- 
The decay per-trigger yield from 7r°'s can then be expressed in terms of the parent 
per-trigger yield. 

2 Constraining Direct Photon-Jet Production 
in Elementary p + p Collisions 

2.1 Isolation Cut Analysis 

In addition to the above statistical subtraction method, standard photon isolation cuts 
were applied event by event in a new analysis in order to dramatically reduce the 
contamination of di-jet events with ir° decay or fragmentation photons in the 7-hadron 
event sample. To be considered isolated, the sum of p^from all tracks and and EMCal 
energies must be < 0.1i? 7 in a cone around the photon of size AR = \J A(j) 2 + Aif = 
0.5. Statistical subtraction of the remaining contribution for isolated ir° production is 
achieved through a modified version of the statistical method above where isolated 7r° 
are used as input to the decay-photon calculation. 

As in the statistical method the analysis is performed as a function of Acf>. We find 
good agreement with the statistical method results as shown in figure HJ Since the sta- 
tistical subtraction method includes, in principle, the contribution from fragmentation 
photon triggers, the agreement between the two methods places a constraint on the 
magnitude of such a contribution. 

2.2 Per trigger yields and Awayside Jet Fragmentation 
Function Analysis 

Results from the isolation cut analysis are shown in figure [21 plotted in terms of the 
fragmentation variables zx = p T 3SOC /Px and xe = %t cos A<p (taking only the compo- 
nent of the associated hadron momentum that is in the same direction as the trigger 
photon). We measure in 6 photon trigger px bins covering the range 5-15 GeV/c. The 
inverse slope parameter of exponential fits to the xe distribution with 7r° and direct 
photon triggers are shown in figure [3l Comparisons with theoretical calculations are 
currently underway. 

2.3 Model Dependent Determination of kx 

Also in pQ, a method for extracting the apparent intrinsic kx itself was described and 
used to extract a value of kx= 2.68 ± 0.35 for 7r°-h correlations. Using isolated direct 
photon triggers, the method is simplified and can be also be used to extract a kx value 
in this channel. Please refer [1] for details of the method. The simplification occurs 
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Figure 1: For two trigger photon bins as indicated on the figure, the awayside per-trigger 
yield vs. A<p for the statistical and isolation cut analysis. The two methods give consistent 
results with the isolation method having much improved precision due to having only a small 
decay photon background subtraction. 




Figure 2: Per trigger yields of direct photon-hadron pairs as a function of Xe/zt for various 
selections of trigger p?. The data have been scaled by factors of 10 for visibility. 
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Figure 3: The inverse slope parameter of exponential fits to the xe distributions for ir° direct 
photon triggers. 



due to the lack of an additional fragmentation variable on the nearside, i.e prt = 
Prt = PT'y- The measured experimental variable p ou t whose distribution for various 
trigger pj> bins is shown in figure 2J is defined as px a si n A</>, or the associated hadron's 
transverse momentum component perpendicular to the trigger photon direction, and 
is proportional to the kx on an event by event basis. 
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Figure 4: p ou4 distributions for various values of trigger px- The data have been scaled by 
factors of 10 for visibility. 

By finding the average p ou t for each trigger photon pt bin as shown in [5] a) one 
can extract the quantity average hr/x where x is just the ratio of the true awayside 
jet momentum to that of the direct photon trigger. This can be extracted from other 
measurements as in [1] and eliminated, leaving a pure measurement of the kj<. For now, 
we simply use the event generator PYTHIA 6.3 to generate direct photon processes 
with initial and final state radiation turned off, and phenomenological hr parameters 
as shown on [5] b) to extract x and make a comparison to PYTHIA, which shows that 
for PYTHIA- like distributions of x a parameter in the vicinity of 2.5-3 GeV gives 
values similar to the data, well-consistent with the value of 2.68 found for di-hadron 
di-jet correlations. 

This implies that the same kx bias effects in the di-hadron correlations exist also in 
the direct photon-jet correlations. This makes the comparison with the perturbative 
QCD calculations, which at only NLO should not contain such kx bias modifications 
even more important, as well as understanding in energy loss models of Au + Au how 




Figure 5: a) y< p 2 out > vs trigger p T for direct photons b) k T /x vs trigger p T for direct 
photons along with calculations from PYTHIA 6.3. 



such a large value of imbalance between the trigger photon and awayside jet might 
alter energy loss interpretations. 

2.4 Direct Measurement of Prompt Photons from Jet Frag- 
mentation 

PHENIX has made further strides towards understanding complications to the LO 
intepretations of direct photon-jet correlations, in studying the contribution of single 
prompt photons that occur in di-jet events, the so-called fragmentation photons pro- 
duced by NLO hard photon radiation or non-perturbative fragmentation processes. A 
direct measurement technique discussed in [10] has been used to measure the angu- 
lar distribution of such photons with respect to trigger hadrons in events where high 
Pt photons are tagged to only be nearby such hard hadrons-an aniz-isolation selec- 
tion. Integrating this distribution one can find the fraction of such photons to total 
hadron-correlated photons from all sources including decay, shown in figure [6] which 
can possibly be related to the fraction of the total direct photon production rate due 
to these fragmentation photons. However due to the restricted kinematic region of the 
measurement, interpretations of this fraction and the angular distribution itself need 
input from sophisticated higher order pQCD calculations (which likely do not yet exist) 
and thus close attention from the theory community. Nonetheless, the measurement 
is exciting because it provides the first step towards making the measurement in the 
Au + Au environment where several interesting predictions of jet-medium enhancement 
of the rate of production of such photons exist. 

3 Constraining Energy Loss theory in Au+Au 
with 7- jet 

In previous conferences |lllll2j PHENIX has presented Au+Au results of the statistical 
subtraction method for direct photon-hadron correlation yields. With the now statis- 
tically improved p + p results (for this section, using the pure statistical subtraction 
method in p+p for comparison, combining Run6 and Run6 statistics) and an expanded 
Pt range for the trigger photons in Au + Au, we can now divide the per-trigger yields 
for Au + Au and p + pin many pt bins making the ratio Iaa = ^direct /^direct f° r an 
awayside integration range of 2ir/5 radians around ir. An example is shown for the 
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Figure 7: Ratio Iaa (see text) for the lowest trigger direct photon p T bin. 

direct photon (trigger) px bin 5-7 GeV/c in figure [7J It is apparent that there are very 
large uncertainties, but most points have positive yield in Au + Au and a value of Iaa 
between and 1 which indicates the basic expectation of suppression of the awayside 
jet. 

We hnd that different sources of uncertainties dominate in different ranges of trigger 
photon and associated hadron hadron px- As either px is increased the statistics 
in the raw correlation functions obviously decrease as the production probabilities 
follow steeply falling spectra. However the combinatoric background from pairs of hard 
particles and soft underlying or otherwise uncorrelated particles is reduced dramatically 
as the px of the associated hadron is increased, having the effect of reducing the 
systematic uncertainty from the underlying event subtraction. Further, as the px of 
the trigger photon is increased there is a larger fraction of direct photon-h pairs in the 
inclusive photon-h pair sample due to the increasing direct photon signal relative to 
the suppressed decay photon background, and thus the higher trigger photon px bins 
also have a reduction in the uncertainty from the subtraction in equation [TJ Because of 
these effects we find that the total uncertainty remains more constant with increasing 
Pt than the loss in statistical precision that the falling production probability would 
normally manifest. For this reason, integrating over large px hadron or trigger px 
bins causes a loss of information since the steeply falling production probability causes 
uncertainties to be dominated by the source corresponding to the low end of the px 
bin. Therefore in order to effectively give the higher px bins more weight, we take the 
plain average of all px,-y,PT,h bins. We call this average the Mean Value Iaa akin to 



1.2 

a 0.8 t 
a 

to 0.6 
> 

c 0.4 

eg 

^ 0.2 



t 1 1 — .— i 1 — | — . 1 — i — | . — i ■ — | — i 1 r- 

Directv-h Run 4 200 GeV Au+Au 0-20% 



PHENIX Preliminary 



2 < PS**" < 10 GeV/c 



10 



a) 



PH ENIX 



12 



14 

T, trigger 



n 
> 

e 

n 



1.6 
1.4 
.3 1.2 
1 

0.8 
0.6 
0.4 
0.2 




: Run 4 200 GeV Au+Au 0-20% #y .„ 
r 5 <p tt _<i5G*v pH - EN(X 

2< p ri <10G«V ■— . * I^Raa 


■ 

pi 






l 




* * t ( 








50 100 1 50 200 250 300 350 

Npar, 



Figure 8: a) Mean value Iaa for each trigger pt bin b) over all px bins for 0-20%, 20-40%, 
and 40-92% centrality bins. 



performing the functional mean value: 

MeanValue Iaa = -r / I A A{PT S ° C )dPT (2) 

Under the further assumption that Iaa remains constant with px this average 
indeed corresponds exactly to the functional mean value, and it is found that this 
assumption is satisfied to what is likely a sufficient degree (considering our large un- 
certainties which would dominate any such uncertainties in the assumption) both in 
measurements of di-hadron correlations [I] and [3] and in most theoretical predictions 
[THE] at sufficiently high pt (associated pr >— GeV/c). 

Results of the Mean Value Iaa are shown in figure El The transition of different 
uncertainty sources from systematic to statistical is even more apparent in figure [8] a) 
for central events, and indeed the Mean Value Iaa seems to favor a constant value 
vs. trigger direct photon pt although the uncertainties are still much too large to rule 
out with any significance rather large possible trigger photon px dependence scenarios. 
Still, again taking the average of all Iaa points together, shown on figure [8] b), now 
also for more peripheral bins, we find that the Mean Value Iaa for the central events 
is significantly positive at a two-sigma level and consistent with the single particle 
suppression level Raa- 

This confirms the basic geometrical picture of energy loss at RHIC |15|.|16j because 
vertices in 7-h correlations, and the trigger direct photon themselves, come from the 
entire collision volume. Therefore the fraction of 7-/1 pairs that are observed without 
significant suppression, come from the same geometrical region that surviving single 
high pt particles do. Dividing this quantity by the unsuppressed direct photon yield is 
analogous to dividing the surviving single particle production rate by the expectation 
from p + p multiplied over the entire production weighted geometric volume of the ini- 
tial distribution of all hard scattering production points, as is done in the construction 
of .R^This means that the away-side suppression in the 7 — jet channel should simply 
reflect the same geometry as the single particle picture and, if geometry plays a domi- 
nant role, give the same suppression level. Since this geometrical picture is believed to 
be a "surface bias" picture where the only di-jets initially produced near the surface 
it could be a confirmation of the surface bias picture, although further comparison 
to the di-hadron Iaa comparing 7-hadron and di-hadron correlations together in the 
same energy loss framework, are necessary to make this statement, in order to rule out 
possible effects in the di-hadron Iaa that are believed not to be from geometry. 

Looking forward to the full release of our Run7 analysis, we also present a first look 
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Figure 9: A0 distribution of direct photon-hadron correlations from Run7 data for 5-7 trigger 
direct photon bin. 



in figure [9] at the lowest trigger px bin using ~ 3.0 billion events, which entices us to 
look across several px bins to look for consistent behavior that may be consistent to 
jet-medium displaced peak as seen already in di-hadron correlations, although at this 
point uncertainties are too large to make a definitive statement with just a single pt 
bin combination. 



4 Conclusions 

Various 2-particle direct photon-hadron correlation yields in p+p and Au+Au collisions 
at ^sjyjy = 200 GeV have been presented. Under the assumption that the suppression 
is nearly px independent using a specific averaging scheme, we find a ratio of Au + Au 
to p + p per-trigger photon yields averaged over all available px bins, Iaa, consistent 
with the single particle suppression level Raa, which can be interpreted as a qualitative 
confirmation of the basic geometrical picture of jet suppression at RHIC. To the extent 
that the prevailing geometric picture is believed to be that of surface bias where only 
jets ejected near the surface dominately contribute to yields, [15], [16], it may be said 
to confirm this picture, however it should be noted that any geometric scenario would 
yield Iaa — Raa- Nonetheless future comparisons of Iaa from 7-jet to that of di-jets, 
may indeed yield some insight into the details of possible surface bias models. 

We have also presented the first event by event isolation cut 2-p correlation results 
at RHIC. The application of the event by event photon isolation cuts in p + p results 
our highest precision measurement yet, and allows for precision studies of the baseline 
fragmentation function D{z), and well as a variable p ou t which is proportional to the 
apparent intrinsic kx, or non-zero transverse momentum of the original collision par- 
tons. With a model dependent extraction method, the average \J < k T > at this y^s 
in p + p is found to be approximately 2.5-3 GeV, consistent with analysis of di-hadron 
(di-jet) correlations pQ. The possible implications of this and also the improved preci- 
sion in the isolated yields warrant detailed comparison with the baseline perturbative 
QCD (pQCD) calculations used in the various models of jet energy loss. Finally, we 
have presented a unique direct measurement of single prompt photons from jet frag- 
mentation, of both angular information and the fraction of these photons in the entire 
photon sample (including decay photons) in the vicinity of the jet cone and for specific 
px cuts. 
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